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THE  NONLINEAR  INTERACTION  OF  INTENSE  LASER 
PULSES  IN  PLASMAS 


I .  Introduction 

The  interaction  of  ultra-high  power  laser  beams^  with  plasmas  is  rich 

2 

in  a  variety  of  wave-particle  phenomena.  These  phenomena  become 

particularly  interesting  and  involved  when  the  laser  power  is  high  enough 

to  cause  the  electron  oscillation  (quiver)  velocity  to  become  highly 

relativistic.  Some  of  the  interesting  laser-plasma  processes  which  are 

3-12 

discussed  include:  a)  relativistic  optical  guiding  of  the  laser  beam, 

b)  the  excitation  of  coherent  radiation  at  harmonics  of  the  fundamental 

13—16 

laser  frequency,  c)  the  generation  of  large  amplitude  plasma  waves 

(wakef ields) ,  d)  frequency  shifts  induced  in  the  laser  pulse  by  plasma 

17  18  19 

waves,  ’  e)  frequency  amplification  using  an  ionization  front  and 

f)  single  particle  acceleration  in  a  laser  pulse. 

In  the  following,  a  fully  nonlinear  ID  model  is  developed  which 

describes  the  self-consistent  interaction  of  intense  laser  pulsps  with 

plasmas.  By  assuming  a  "quasi-static"  cold  fluid  plasma  response,  a  set  of 

coupled  nonlinear  equations  is  derived  for  the  vector  potential  of  the 

radiation  field  and  for  the  electrostatic  potential  of  the  plasma.  The 

quasi-static  approximation  assumes  that  in  a  frame  moving  at  the  speed  of 

light,  the  plasma  fluid  sees  a  nearly  steady-state  radiation  field  (after 

the  transients  have  decayed  away).  The  resulting  nonlinear  equations  are 

used  to  examine  various  laser-plasma  interaction  phenomena.  The  important 

2  3  20 

issue  of  laser-plasma  instabilities  ’  ’  is  not  addressed  in  this  paper. 
Instabilities  will  certainly  limit  the  laser  ]~“?-''gation  distance  and, 
therefore,  will  have  a  profound  effect  on  some  he  potential 
applications  discussed  in  this  paper. 

Manuscript  approved  September  20.  1989. 
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Relativistic  optical  guiding  is  a  result  of  the  relativistic 

quiver  motion  of  the  electrons  by  the  laser  field.  Analysis*^  has  shown 

that  as  the  laser  power  exceeds  a  critical  threshold,  diffraction  can  be 

overcome,  resulting  in  optical  guiding  of  the  laser  pulse.  Previous 

analysis  of  relativistic  guiding  have  included,  for  the  most  part,  only 

3-12 

the  transverse  electron  motion  in  the  plasma  response  current. 

Relativistic  guiding  was  believed  to  occur  on  a  fast  time  scale  (on  order 
of  the  inverse  laser  frequency).  However,  the  present  analysis  finds  this 
not  to  be  the  case.  In  the  following  nonlinear  analysis  of  relativistic 
optical  guiding,  the  electron  density  response  and  longitudinal  electron 
motion  are  included  self-consistently .  It  is  shown  that  for  short  laser 
pulses  (pulse  lengths  less  than  a  plasma  wavelength),  the  combined  effects 
of  the  plasma  density  response  and  the  longitudinal  motion  significantly 
reduce  the  relativistic  guiding  effect.  It  is  found  that  relativistic 
guiding  occurs  only  for  long  pulses  with  slow  rise  times  (greater  than  an 
inverse  plasma  frequency). 

As  the  quiver  motion  of  the  electrons  in  a  linearly  polarized  laser 

field  becomes  highly  relativistic,  the  plasma  response  current  will  develop 

harmonic  components.  The  harmonic  content  of  the  response  current  density 

can  lead  to  the  excitation  of  coherent  radiation  at  harmonics  of  the 

fundamental  laser  frequency.  It  is  shown  that  harmonic  generation  is  more 

effective  for  short  laser  pulses  (pulse  lengths  less  than  a  plasma 

wavelength)  than  for  long  pulses  with  slow  rise  times. 

An  intense,  short  pulse  laser  (pulse  lengths  near  the  plasma 

wavelength)  interacting  with  a  plasma  can  generate  large  amplitude  plasma 

wave  wakefields.  The  excited  Wakefields  may  be  used  to:  i)  accelerate  a 

13—16 

trailing  electron  bunch  (laser  Wakefield  acceleration),  ii)  optically 
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guide  a  trailing  laser  pulse  and  iii)  enhance  the  coherent  harmonic 

radiation  generated  by  a  trailing  laser  pulse. 

Three  other  laser-plasma  interaction  phenomena  are  briefly  discussed. 

The  first  describes  how  frequency  shifts  are  induced  in  a  laser  pulse  by  a 
17  18 

plasma  wave,  ’  the  second  discusses  how  an  ionization  front  may  be  used 

19 

to  amplify  the  frequency  of  a  laser  pulse  and  the  third  describes  how  a 
single  particle  may  be  accelerated  by  the  passage  of  a  laser  pulse. 


II .  Nonlinear  Formulation 

The  ID  fields  associated  with  the  laser-plasma  interaction  can  be 

described  by  the  transverse  vector  and  scalar  potentials,  A^(z,t)  and 

#(z,t),  respectively  (see  Fig.  1).  In  what  follows,  we  use  the  Coulomb 

gauge,  i.e.,  V-A  =  0  which  implies  A  =  0.  The  transverse  polarization  of 
**  z 

the  laser  field  is  arbitrary.  The  normalized  potentials  satisfy  the 
following  equations, 


y 4  ■  kp  k  *  -  kp  k  - /Y ’ 


(la) 


a2* 


3z 


2-t 


(lb) 


2  2 

where  a(z,t)  =  |e|A^/moc  ,  <fr(z,t)  =  |e|$/moc  ,  kp  =  «p/c,  <0p  = 

2  1/2 

( 4 n | e |  no/mQ)  is  the  ambient  plasma  frequency,  n(z,t)  is  the  plasma 

density,  nQ  is  the  ambient  plasma  density,  =  v^/c  =  a/y  is  the 

2  2  -1/2 

normalized  transverse  plasma  fluid  velocity  and  y  =  (1  -  0  -  g^)  =  (1 

2  1/2  2  1/2 

+  a  )  /(I  -  0Z)  is  the  relativistic  mass  factor.  In  obtaining  the 

right-hand  side  of  (la)  we  used  the  fact  that  the  transverse  canonical 
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momentum  is  invariant  and  prior  to  the  laser  pulse  interaction  the  plasma 
is  assumed  stationary. 

The  fluid  quantities  n,  3z  and  y  are  assumed  to  satisfy  the  cold 
relativistic  fluid  equations  which  can  be  written  in  the  form, 

If  ♦  c  h  <“V  -  °-  <2a> 


d t  --  - <1/2)1'-2  (c  h  *  ^  It)3-2  +  T  <>  -  &  ¥z'  <2b> 


If  ■  c6z  H  *  It  a~2/3t’  <2c> 

2 

where  d/dt  =  3/3t  +  c3z3/3z.  The  term  containing  a  in  Eqs.  (2b)  contains 
the  ponderomotive  force.  Thermal  effects  may  be  neglected  provided  the 
electron  quiver  velocity  is  much  greater  than  the  electron  thermal 
velocity.  Also,  the  ions  are  assumed  to  be  stationary. 

It  proves  convenient  to  perform  an  algebraic  transformation  from  the 
laboratory  frame  independent  space  and  time  variables  (z,t)  to  the 
independent  variables  (?.,t),  where  £  =  z-c3tt  and  x  =  t.  Here  Pt  = 
vt/c  =  1  is  the  normalized  transformation  velocity.  To  transform  Eqs.  (1) 
and  (2)  from  z,t  to  ?,,x  variables  we  note  that  3/3z  =  3/3£,  and  3/3t  =  3/3x 
-  c0t3/3f,.  Using  these  transformations,  Eqs.  (1)  and  (2)  become, 

,  9  a2  2(3  .2  9  ,2  .  k?- 

k  ^2)5 '  y  r3-  (33> 


q.  k2  1), 

8t2  »  1 


It  [n(6r  v]  '  c  H' 
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(3d) 


h  [r  (1-w  -  *]  -  -  £1? 

22  21/2 
where  y  =  l/(l-0t).  Equations  (3a-d),  together  with  y  =  (1  +  a  )  / 

2  1/2 

(1  -  3z)  ,  form  a  complete  set  of  fully  nonlinear,  relativistic,  cold 

fluid  equations  which  describe  the  ID  laser-plasma  interaction.  The  ID 

model  is  valid  as  long  as  the  radiation  spot  size  is  large  compared  to  the 

plasma  wavelength,  i.e.,  r  >>  X  =  2n/k  . 

s  p  p 

In  the  following  analysis,  we  choose  v  =  c,  since  v  =  c  lies  between 

the  group  velocity  and  phase  velocity  of  the  laser  pulse.  Furthermore,  we 

take  the  laser  pulse  to  propagate  in  the  positive  z  direction  and  to  cross 

the  z  =  0  plane  at  t  =  0.  The  £  =  0  plane  defines  the  front  of  the  pulse 

which  extends  into  the  £  <  0  region.  In  the  £,T  frame  we  are,  therefore, 

interested  only  in  the  region  where  £  <  0  since  for  (  >  0  we  have  a  =  0, 

n  =  n  ,  0  =0  and  y  =  1. 

o  z 


a)  Quasi-Static  Approximation 

Equations  (3)  can  be  greatly  simplified  by  noting  that  in  the  speed  of 

light  frame  (0=1)  we  expect  that  under  certain  conditions  a  quasi-static 

state  will  exist  in  the  macroscopic  plasma  quantities,  n,  0z  and  y.  That 

is,  if  the  laser  pulse  is  sufficiently  short,  the  fields  a  and  4>  which 

drive  the  plasma  are  expected  to  change  little  during  a  transit  time  of  the 

plasma  through  the  laser  pulse.  From  (3a)  we  find  that  the  envelope  of  a 

changes  on  a  characteristic  time  T  ~  2y|n  /n|(co/«  )/«  ,  where  w  is  the 
°  e  1  o  1 '  p'  p' 

laser  frequency.  Since  we  will  henceforth  assume  that  «  >>  the 
radiation  envelope  changes  on  a  time  scale  which  is  long  compared  to  a 
plasma  period.  If  the  laser  pulse  duration,  t^,  is  small  compared  to  t  , 
i.e.  T,  <<  T  ,  then  the  quasi-static  approximation  is  valid.  In  addition, 

Li  6 
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the  validity  of  the  ID  model  requires  that  the  laser  beam  diffraction  time 

2 

(transverse  spreading  time),  =  nrs/(Xc),  be  long  compared  to  Tg.  This 

condition  is  satisfied  as  long  as  r  >>  X  .  More  formally,  the  quasi- 

s  p 

static  approximation  (neglecting  the  right-hand  side  of  (3c, d))  implies 


1  a_ 

c  3t  , 


o 

ndE,' 


<<  n  , 
o 


(4a) 


«  1. 


In  this  case,  Eqs.  (3c, d)  can  be  integrated  to  give 


(4b) 


n(l-3  )  =  n  , 
v  Kz '  o’ 


(5a) 


Y(l-Pz)  -  ♦  =  1. 


(5b) 


Using  (5a, b)  together  with  the  expression  for  y>  the  coupled  field 
equations  become 


(l  i_  c-2  l_)  _  k2 

lc  as,  a-cj  3t  "  p  uy 


(6a) 


^-4  +  k2(l  +  <j>)  2 ( 1  +  4*/ 2 ) <J>  =  k2( l+4»)  2  a2/2, 


(6b) 


ar 


where  the  plasma  quantities  in  terms  of  the  fields  are 

n/nQ  =1+|  ((l+a2)/(l+*)2  -  l),  (7a) 

Y  =  (l+a2  +  (U*)2)/ [2(1+*)) ,  (7b) 

3z  =  (l+a2  -  (l  +  *)2)/(l+a2  +  (1+*)2).  (7c) 
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The  above  expressions  for  the  fields  and  plasma  quantities  are  fully 
nonlinear  and  hold  for  arbitrary  polarization  of  the  laser  pulse. 
Consistent  with  the  quasi-static  assumption,  a  number  of  important  points 
can  be  made  concerning  intense  laser  pulse  propagation  in  plasmas. 


Ill .  Laser-Plasma  Interaction  Phenomena 
a)  Relativistic  Optical  Guiding 

The  nonlinear  index  of  refraction  of  the  laser  beam  within  the  plasma 
determines,  among  other  things,  the  optical  guiding  properties  of  the 
plasma.  For  the  purpose  of  the  present  discussion,  we  assume  that  the 
laser  field,  a,  is  given  by 

a  =  aLa,x)elkC/2  +  c.c.,  (8) 


where  a L  represents  the  complex  amplitude  and  k  is  the  wavenumber.  The 
characteristic  spatial  variation  in  the  laser  envelope,  |aL|,  is  assumed  to 
be  of  order  L  and  is  long  compared  to  the  laser  wavelength,  X,  i.e., 
3|aL|/3£  =  |aJ/L  «  k|aj  =  2n|aL|/A. 

Using  the  representation  in  Eq.  (8)  we  find  from  Eq.  (6a)  that  the 
refractive  index,  h  =  ck/co,  is 

=  1  -  \  (kp/k)2/(l  +  *s),  (9) 


where  <|>  is  the  slow  part  of  the  scalar  potential  and  k  <<  k.  In 
s  p 

obtaining  Eq.  (9),  we  replaced  (l/(l+$))g  with  l/(l  +  <J>g),  which  is  valid  as 

long  as  kp  <<  k  (this  implies  that  |^|  <<  |<f>s|,  where  <|>^  is  the  rapidly 

varying  part  of  4>) .  There  are  two  cases  that  can  be  considered,  depending 

on  the  envelope  scale  length,  L,  compared  to  the  plasma  wavelength,  X  . 

P 
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i)  Long  Rise  Time  Pulse 

In  this  limit,  L  >>  X  =  2nc/u>  ,  the  first  term  on  the  left-hand  side 

P  P 

2  1/2 

of  (6b)  can  be  neglected  and  <p  can  be  approximated  by  <$>  =  ( 1  +  | a.  |  / 2 )  -1. 

S  S  Li 

The  refractive  index  becomes 

h  =  1  -  \  (X/Xp)2(l+|aL|2/2)~1/2.  (10) 

Although  the  present  analysis  is  ID,  we  expect  that  for  a  slowly 
varying  transverse  laser  profile  the  index  of  refraction  will  depend  on  the 
transverse  coordinates  through  the  laser  amplitude  |a^|.  Since  the  actual 
transversely-bounded  laser  beam  amplitude  falls  off  transversely,  i.e., 

9  la.  l/9r  <  0,  so  will  the  refractive  index,  i.e.,  9h|/3r  <  0.  The  negative 

Li 

transverse  gradient  of  the  refractive  index  can  lead  to  optical  guiding. 

Since  the  refractive  index  is  a  function  of  the  laser  amplitude,  the 
condition  for  optical  guiding  places  a  lower  limit  on  |a^|.  It  is  well 
known  that  if  the  refractive  index  is  of  the  form  given  by  Eq.  (10),  a 
critical  laser  power  necessary  for  relativistic  optical  guiding  exists  and 
is  given  by6~9  Pcrit  =  17.4(Xp/X)2[GW] . 

ii)  Short  Pulse 

Consider  the  short  laser  pulse  limit,  L  <  Xp.  When  |<|>|  <<  1,  Eq.  (6b) 
can  be  solved  for  an  arbitrary  laser  field  a(C)> 

k  ° 

4>  =  -f  J  a2(£')sin(kp(*;'-S))d*;',  (11) 

where  the  boundary  conditions  <f>(£=0)  -  3<KS=0)/3£  =  0  have  been  used. 

If  the  pulse  envelope  is  given  by  a^  =  a^Qsin( ix£,/L)  for  -L  <  £  <  0  and 
a.  =  0  otherwise,  we  find  that  the  scalar  potential  within  the  laser  pulse 

Li 

is  given  by 
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4>  =  (aLo/8){1_(kp  '  4ti2/l2)  JkpCOs(2n£/L)  _  (4n2/L2)cos(kp£.)]|,  (12) 

2 

where  terms  of  Older  (X/Xp)  <<  1  have  been  neglected.  For  L  <<  Xp,  this 
gives  $  =  (a.  k  /4)2  g(0>  where  g(£)  =  £2-2(L/2it)2[  1-cos(2ji£/L)  1 .  Notice 

S  LO  p 

2 

that  for  L  <<  X  ,  $  is  maximum  at  E,  =  L  where  <J>  =  (a.  k  L/4)  .  Also 

p  s  s  '  Lo  p 

notice  that  even  for  |a,  |  >  1,  the  assumption  that  <t>  <<  1  is  still  valid 

LO  S 

as  long  as  L  <<  X  .  The  index  of  refraction  in  the  short  pulse  limit  is, 
therefore, 


1  -  \  ( X/ X  ) 2 / [  1  i-  (it/2)2(aLo/X  )2g(0], 


(13) 


where  -L  <  £  <  0. 

In  the  short  pulse  limit,  the  fact  that  <t>s  <<  1  implies  that  the 

optical  guiding  effect  is  reduced  significantly,  by  more  than  the  factor 
o  2 

(JlV2)(L/X  )  <<  1.  The  critical  power,  therefore,  is  increased  by  the 

inverse  of  this  factor  and,  in  addition,  the  degree  of  guiding  varies  along 
the  pulse.  Hence,  it  is  unlikely  that  relativistic  optical  guiding  can  be 
effectively  utilized  in  short,  L  <  Xp,  laser  pulses. 

Although  it  may  appear  that  a  long  laser  pulse  may  undergo  guiding, 
assuming  the  various  laser-plasma  instabilities  can  be  controlled,  the 
front  of  the  pulse  will  diffract.  Initially,  that  portion  of  the  head  of  a 
long  rise  time  pulse  in  which  the  local  power  is  less  than  will 

diffract.  Once  this  portion  has  diffracted  away,  the  pulse  will  exhibit 
"short  pulse"  diffractive  behavior,  i.e.,  the  front  region  (~  Xp)  will 
continue  to  diffract.  The  erosion  of  the  front  of  the  pulse  due  to 
diffraction  will  propagate  back  through  the  body  of  the  pulse.  The  erosion 
velocity  back  through  the  body  of  the  pulse  (in  the  £  =  z-ct  frame)  may  be 
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estimated  by  v£  =  (Xp/zR)c,  where  zR  =  itrVX  is  the  vacuum  Rayleigh  length, 
b)  Harmonic  Excitation 

The  nonlinearities  associated  with  the  plasma  waves  can  provide  a 
source  for  the  generation  of  coherent  radiation  at  harmonics  of  the  laser 
frequency.  To  examine  this  process  we  generalize  Eq.  (8)  and  represent  the 
full  radiation  field  by 


=  a.  ( £,  't)e^*c^’/2  +  c .  c . , 


(14) 


where  i  =  1,2,3,...  and  a^  =  aL(£)  is  the  envelope  of  the  dominant 
fundamental  laser  pulse,  |aL|  >>  ja^|,  for  t  >  2. 

It  is  clear  from  the  right-hand  side  of  (6a)  that  harmonic  excitation 
is  solely  due  to  the  excitation  of  the  fast  part  of  <j>.  In  particular, 
since  the  fundamental  component  of  the  radiation  field  dominates,  the  fast 
part  of  the  scalar  potential  is 

4>f  =  -  (kp/4k)2a£(l  +  <y“2cos(2kO,  (15) 

where  we  have  taken  the  fundamental  laser  pulse  to  be  aL(E,)cosk£.  The 
source  term  in  Eq.  (6a)  becomes 

S  =  k2aL(l  +  4»s+<j>f)“1cos(kO  =  k2aL(l  +  <f>s)"1cos(kOE(«cos(2kO)m,(16) 
where  <f>f  is  given  by  (15),  a  =  ((k  /4k)a.  )2(l+$  )-2  and  m  =  0,1,2,... 

L  p  1_»  S 

As  an  illustration  we  consider  the  excitation  of  third  harmonic 
radiation  (3«).  Since  |4>^|  <<  |<|>s|,  the  third  harmonic  component  of  the 
source  is 

S3  =  (k2/2)(kp/4k)2aJ(l  +  4»s)~4cos(3kJ;),  (17) 
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where  we  have  used  the  m  =  1  term  in  Eq.  (16).  Substituting  (17)  into  the 
right-hand  side  of  Eq.  (6a)  and  solving  for  the  third  harmonic  field  we 
find  that 


la3 1 =  j  4aLy>  <18> 

where  x  is  the  laser-plasma  interaction  time.  The  ratio  of  the  third 
harmonic  power  to  the  fundamental  laser  power  is 


2 

P3/Pl  '  [<V4X)3(1.*s)-4^.t]  , 


(19) 


10  2 

where  P.  =  2.15  x  10  (r  a,/X)  [Watts]  (for  a  Gaussian  transverse  laser 

1  S  L» 

profile) . 

Equation  (16)  shows  that  the  generation  of  harmonics  is  a  strong 

function  of  the  plasma  wakefield,  described  by  $  ,  in  the  region  of  the 

s 

fundamental  laser  pulse.  For  a  single  long  pulse,  large  amplitude  laser, 

|aL|  >>  1,  the  slow  part  of  the  scalar  potential  is  <f>s  ~  |aL|//2.  In  this 

case,  the  harmonic  content  of  the  source  term  in  Eq.  (16)  is  exceedingly 

3  2  2 

small,  and  for  the  third  harmonic,  we  have,  P^/P-^  ~  (MXMX^)  c*>pX/  | aL  |  )  . 

2 

Taking  x  to  be  a  diffraction  time,  x^  =  itrs/Xc  and  rg  >  X^,  the  third 

A  A  i  i  A 

harmonic  power  becomes  P^/P^  ~  (rt/2)  (X/X^)  /(4|a^|  ).  In  the  case  of  a 
short,  L  «  Xp,  large  amplitude  laser  pulse,  <|>s  <  |aL|  (HL/2X  )  ,  the  third 
harmonic  power  is  P^/P j  =  (  |a^  |^(X/4Xp)^(0pX)^  =  |aL  |^( ix/4)^(X/Xp)^/4  for  x 
=  x^  and  rg  >  Xp.  For  |aL|  >>  1,  a  short  pulse  is  more  efficient  than  a 
long  pulse  for  harmonic  generation. 

Harmonic  generation  can  be  enhanced  in  regions  of  a  plasma  wave  where 
4>s  <  0.  This  can  be  achieved  by  positioning  a  short  laser  pulse  at  the 
appropriate  location  in  a  large  amplitude  plasma  wave  having  a  phase 
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velocity  close  to  the  speed  of  light.  The  large  amplitude  plasma  wave  can 

13—16 

be  the  wakefield  generated  by  a  laser  pulse  or  by  an  electron  beam 

i  21 

pulse. 


c)  Wakefield  Generation 

The  self-consistent  evolution  of  the  nonlinear  plasma  wave  can  be 
studied  by  numerically  solving  Eqs.  (6a, b).  Figures  2a, b  show  the  plasma 
density  variation  8n/nQ  =  n/nQ  -  1  and  the  corresponding  axial  electric 
field  Ez  for  a  laser  pulse  envelope  given  by  a^  =  a^Qsin(  rtJ,/L)  for 
-L  <  J,  <  0.  In  this  figure,  L  =  X  =  0.03  cm,  X  =  10  p,  and  a)  a,  =  0.5 

P  LO 

and  b)  a^Q  =  2.  The  steepening  of  the  electric  field  and  the  increase  in 

22-24 

the  period  of  the  wakefield  are  apparent  for  the  highly  nonlinear 

situation  shown  in  Fig.  2b  ( |aLo  |  =  2)  as  compared  to  the  slightly 
nonlinear  case  shown  in  Fig.  2a  (  |aLo  |  =  0.5).  Figure  3  shows  that  the 
electrostatic  potential  <f>  is  predominantly  slowly  varying  within  the  laser 
pulse  even  though  8n/no  has  rapidly  varying  components.  In  addition,  the 
figure  shows  that  <f>  can  be  negative  behind  the  laser  pulse. 

Qualitative  aspects  of  nonlinear  wakefield  generation  may  be  examined 


by  considering  a  circularly  polarized  laser  pulse  with  a  square  pulse 
profile,  i.e.,  a^  =  aLo  for  -L  <  £  <  0  and  a^  =  0  otherwise  [also,  see 
Ref.  25].  For  this  case,  Eq.  (6b)  may  be  solved  analytically  in  terms  of 
elliptic  integrals.  In  particular,  one  can  show  that  the  optimal  pulse 


length  L  =  L  (which  corresponds  to  the  maximum  wakefield  amplitude)  is 

2 

given  by  L  =  2y^QE(p)/kp  -*  2y^Q/kp  for  y^Q  >>  1,  where  E(p)  is  the 

2  2  2  2 

complete  elliptic  integral  of  the  second  kind,  p  =  (y^Q-l )/y^0  and  y^Q  = 

2  2  2 
l+a^Q.  This  gives  a  wakefield  behind  the  pulse  for  which  y^Q  >  1  +  4>S  >  l/y^Q 


and  the  maximum  axial  electric  field  of 

( y?  -l)/y,  ,  where  E  =  I e  I E  /(m  c^k  ). 
lo  '  lo  1  1  z  o  p 


this  wakefield  is  given  by  E 

'  max 

The  nonlinear  wavelength  of  the 
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NL  2  2 

wakefield  is  given  by  X  =  4y,  E(p  )/k  -»  4y,  /k  for  y,  >>  1,  where  p  = 

p  lo  o  p  lo  p  lo  o 

The  large  amplitude  axial  electric  fields  associated  with  the  plasma 

waves  can  be  utilized  to  accelerate  an  injected  beam  of  electrons  to  high 
13—16 

energies  (LWFA).  In  the  region  where  Sn/nQ  <  0,  the  transverse 

profile  of  the  plasma  wakefield  can  lead  to  a  negative  transverse  gradient 

26 

of  the  refractive  index,  3h/8r  <  0.  Equation  (9)  indicates  that  a 
properly  phased  trailing  laser  pulse  (located  at  a  maximum  in  <t>s)  may, 
therefore,  be  optically  guided.  From  Eq.  (16),  it  is  seen  that  harmonic 
generation  can  be  substantially  enhanced  in  a  properly  phased  (located  at  a 
minimum  of  <f>s)  short  trailing  laser  pulse  propagating  in  the  wakefield 


generated  by  a  leading  laser  pulse  when  -1  <  <J>  <  0.  In  addition,  the  sharp 

axial  gradient  in  &n/nQ  for  a  highly  nonlinear  plasma  wakefield  could 

17  18 

induce  large  frequency  shifts  in  a  short  trailing  laser  pulse.  ’ 


d)  Laser  Pulse  Frequency  Variations 

As  the  laser  pulse  propagates  in  the  plasma,  the  excited  plasma  wave 

17  18 

can  modify  the  laser  frequency.  '  To  analyze  this  effect  we  use  the 
representation  in  Eq.  (8)  and  set  aL(E.,t)  =  |a^( £)  |exp( iij>( £,  t) ) ,  where 
|a^(&)|  denotes  the  laser  pulse  envelope,  and  4»(£,t)  is  real  and  denotes 
the  slowly  varying  phase  shift.  Substituting  (8)  into  (6a)  we  find  that 


T 

J  (1  -  4>/(l  +  *>)dT'. 


o 


(20) 


The  laser  frequency  is  w^S.t)  =  -9(k£  +  vj»)/3t  =  ck  +  (c3/3£,  -  3/3t)<j/  = 

2 

ck(l  +  (kp/k)  /2)  +  Aw,  where  Aw  is  the  frequency  shift  due  to  the  plasma 
wave, 
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Aw  = 


(21) 


T 


r  (kP/k)2(u*  ~ c  h  I  u* dT') 


0 

For  completeness,  we  note  that  the  wavenumber  is  given  by  k^(£,  t)  = 
3(k£  +  ^)/3£,  =  k  +  Ak,  where 


X 

- 1  <yk)2  c  ft  J  iVt' • 

o 


(22) 


and  the  refractive  index,  t|  =  ck^/w^,  is  given  by  Eq.  (9).  The  validity  of 
(21)  is  based  on  the  assumption  that  is  slowly  varying;  this  breaks  down 
when  the  frequency  shift  approaches  the  laser  frequency,  i.e.,  |Aw|  ~  ck. 

For  |<f>|  <<  1,  <f>  is  given  by  Eq.  (12)  and  the  laser  frequency  shift  for 
WpT  >>  1  is  given  by 

o 

Aw  =  -  ck(kp/k)2  (wpx/4)kp  |  |aLa')|2coskpa'  -  Od£'.  (23) 


For  a  short  pulse  laser  (L  <<  Xp),  assuming  an  interaction  time  x  =  x^,  and 

spot  size  r  >  X  ,  the  frequency  variation  within  the  pulse  is  given  by 
s  -  p 

3  2  2 

Aw/ w  =  (n/2)  (X£/Xp)a^o  <<  1.  Frequency  variations  induced  on  a  short 
pulse  laser  propagating  in  a  plasma  wave  generated  by  a  leading  laser  pulse 
can  also  be  analyzed.  ^ 


e)  Frequency  Amplification  Using  an  Ionization  Front 

A  laser-induced  ionization  front  may  also  be  used  to  upshift  the 

19 

frequency  of  a  laser  pulse.  An  intense  pump  laser  with  a  short  rise  time 
(<  1.0  psec)  may  be  used  to  ionize  a  gas  as  it  propagates.  This  creates  a 
moving  ionization  front  in  which  the  plasma  density  goes  from  zero  to  large 
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21  -3 

values  (-  10  cm  )  within  a  very  short  distance  (<  0.01  cm).  This  sharp 

plasma  density  gradient  propagates  at  the  group  velocity  of  the  pump  laser, 

2  2  1/2 

which  may  be  approximated  by  v  /c  =  (1  -  w  /«  )  ,  where  «  is  the 

go  p  0  o 

frequency  of  the  pump  (ionizing)  laser  and  eo^  is  the  plasma  frequency  of 

the  ionized  gas.  One  possible  method  for  upshifting  the  frequency  of  a 

laser  pulse  is  the  following.  A  secondary  short  (<  0.1  psec)  laser  pulse 

propagates  nearly  parallel  with  the  pump  laser.  By  phasing  the  secondary 

pulse  such  that  it  "rides"  the  large  density  gradient  of  the  ionization 

front,  the  frequency  of  the  secondary  pulse  may  be  upshifted  in  much  the 

17  18 

same  manner  in  which  a  plasma  wave  may  be  used  '  to  upshift  the 

frequency  of  a  laser  pulse.  For  a  linear  gradient  in  the  ionization  front 

moving  at  velocity  c,  it  can  be  shown  that  the  frequency  of  the  secondary 
laser  pulse  is  given  by 


where  w(0)  is  the  initial  frequency  of  the  secondary  laser  pulse, 

3(n/nQ)/3C  =  -  1/d,  d  is  the  spatial  gradient  length  of  the  ionization 

front  and  z  is  the  interaction  length.  As  an  example,  consider  a  KrF  laser 

21  -3 

(X(0)  =  0.26  pm),  an  ionization  density  of  nQ  =  10  cm  (X^  =  1.0  mm) 

with  d  =  0.01  cm  and  z  =  1.0  m.  This  gives  co(z=l  m)  =  26  w(0)  or  X  =  0.01 

pm.  Notice  that  it  may  be  possible  to  tune  the  upshifted  frequency  by 

varying  the  interaction  length.  However,  due  to  the  difference  between  the 

laser  pulse  velocity  and  the  ionization  front  velocity,  the  laser  pulse 

will  eventually  slip  out  of  the  region  of  the  ionization  front,  thus 

27 

limiting  the  interaction  distance. 

19 

Alternatively,  more  dramatic  "amplifications"  in  the  frequency  may 
be  achieved  by  considering  a  secondary  laser  pulse  propagating  nearly 
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collinear  to  the  pump  laser,  but  in  the  opposite  direction.  In  this  case 

the  ionization  front  of  the  pump  laser  appears  as  a  relativistic  plasma 

mirror  moving  towards  the  secondary  pulse.  Hence,  part  of  the  incident 

secondary  pulse  will  reflect  off  this  ionization  front.  The  frequency  of 

the  reflected  radiation,  cog,  will  be  relativistically  doppler  shifted, 

2 

w  =  4(»  /ft)  )  ft).,  where  ft),  is  the  frequency  of  the  incident  radiation, 
s  o  p  1  X 

This  frequency  amplification  may  be  quite  large,  and  one  can  envision  «  to 

s 

be  in  the  x-ray  regime.  The  power  of  the  reflected  radiation,  however, 

will  most  likely  be  small  for  two  reasons:  i)  The  power  of  the  incident 

radiation  needs  to  be  sufficiently  low  so  as  not  to  appreciably  ionize  the 

gas,  and  ii)  the  ionization  front  may  be  a  "poor"  mirror,  the  reflectivity 

of  which  is  dependent  on  the  profile  of  the  plasma  density  as  well  as  on 

the  ratio  ft). /ft)  .  It  should  be  noted  that  an  electron  beam  may  also  be  used 
1  p 

in  place  of  the  pump  laser  to  create  an  ionization  front. 

f)  Single  Particle  Acceleration 

As  a  final  point  we  note  that  electrons,  initially  at  rest,  can  be 

accelerated  axially  by  experiencing  only  the  transverse  laser  fields.  The 

axial  acceleration  is  due  to  the  v  x  B  force  on  the  electrons.  In  the 

absence  of  a  plasma,  the  scalar  potential  vanishes,  <(>  =  0.  Equation  (7b) 

shows  that  after  experiencing  the  fields  of  a  laser  pulse,  an  electron, 

2  2 

initially  at  rest,  will  acquire  a  final  energy  given  by  =  (a^/2)mQc  , 
where  a^  is  the  final  value  of  the  laser's  vector  potential.  This  is 
reasonable  since  a^  is  proportional  to  the  area  under  the  electric  field. 
From  Eq.  (7c)  and  (7d)  we  find  that  the  final  ratio  of  the  magnitude  of  the 
axial  to  transverse  velocity  is  =  (2+a^)  /(2(l+a^)  |a^  |)  -*  a^/2, 

where  the  second  expression  is  valid  for  |a^  |  >>  1. 
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IV.  Discussion  and  Conclusions 


Based  on  a  ID  nonlinear  quasi-static  model,  a  number  of  different 
laser-plasma  interaction  phenomena  have  been  discussed.  Relativistic 
optical  guiding  is  shown  to  depend  strongly  on  the  laser  pulse  duration. 

In  the  long  pulse  regime,  optical  guiding  requires  a  minimum  level  of  total 
laser  power,  i.e.,  the  critical  power.  However,  the  leading  portion  of  the 
pulse  will  experience  diffraction.  In  the  short  pulse  regime,  relativistic 
guiding  effects  are  greatly  diminished  by  the  density  response  and  the 
longitudinal  motion  of  the  electrons.  As  the  electron  quiver  motion 
becomes  highly  relativistic,  the  plasma  response  current  develops 
harmonics.  This  process  is  analyzed  and  it  is  shown  that  coherent  harmonic 
generation  is  more  effective  for  short  pulses  than  it  is  for  long  pulses 
with  slow  rise  times.  The  generation  of  nonlinear  plasma  wakefields  by 
intense,  short  pulses  was  examined,  and  wave  steepening  is  observed  for 
very  intense  laser  pulses.  Various  applications  of  the  plasma  wakefields 
may  be  possible,  including  i)  the  acceleration  of  a  trailing  electron  bunch 
(laser  wakefield  acceleration),  ii)  optical  guiding  of  a  trailing  laser 
pulse  and  iii)  enhancing  the  coherent  harmonic  radiation  generated  by  a 
trailing  laser  pulse.  In  addition,  the  ID  nonlinear  model  has  been  used  to 
calculate  the  frequency  shifts  induced  in  a  laser  pulse  by  a  plasma  wave, 
as  well  as  single  particle  acceleration  by  a  laser  pulse.  The  possibility 
of  frequency  amplification  by  using  a  relativistic  ionization  front  may 
lead  to  a  source  of  coherent  x-rays.  The  viability  of  this  process  will  be 
a  subject  of  future  research. 
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the  pulse  is  at  £  =  0.  In  this  frame  the  plasma  flows  from  right 


DISTANCE  £ 

Density  variation  Sn/nQ  =  n/nQ  -  1  and  axial  electric  field  in 
GeV/m  for  a  laser  pulse  located  within  the  region  -L  <  E,  <  0, 
where  L  =  X  =  0.03  cm,  and  a)  a,  =  0.3  and  b)  a.  = 


2.0. 


(GeV/m) 
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Fig.  2  (Con't)  Density  variation  6n/no  =  n/nQ  _  1  and  axial  electric  field  Ez  in 
GeV/m  for  a  laser  pulse  located  within  the  region  -L  <  £,  <  0, 
where  L  =  =  0.03  cm,  and  a)  a^0  =  0.5  and  b)  a^Q  =  2.0. 
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Fig.  3.  Electrostatic  potential  <t>  for  a  laser  pulse  located  within  the 
region  -L  <  £  <  0,  where  L  =  =  0.03  cm  and  aLo  =  2.0. 
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